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We show that in the multi-Higgs extension of minimal supersymmetric standard model, which has 
N-pairs of Higgs doublets (called NHMSSM), it is impossible to break CP spontaneously, if we do 
not allow for fine-tuning relations between parameters. The result holds true even in the presence 
of spontaneous R-parity breaking. 
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Evidence for the existence of fundamental forces that distinguish between matter and antimatter (CP violating 
forces) is abundant, both in the laboratory and in cosmology. In the laboratory, there are the celebrated discoveries of 
CP violation in K- and B-systems whereas in cosmology, CP violation is an essential ingredient in our understanding 
of a fundamental mystery of the Universe, i.e., the asymmetry between its matter and antimatter content. While 
this information has been accumulating for over half a century, an understanding of the origin and nature of these 
forces responsible for them has eluded physicists. What is known is that in the standard model, CP violation can 
be parameterized in terms of a single phase in the Cabibbo-Kobayashi-Maskawa quark rotation matrix. However, we 
O ■ do not know where the phase comes from. For instance it is not known whether it resides as an intrinsic phase in 
the interactions of Higgs bosons responsible for particle masses (Yukawa interactions) or it arises dynamically in the 
qj- ground state of the theory, even though all interactions in the model are CP symmetric prior to symmetry breaking. 
The latter class of models go by the name of "spontaneous CP violation" (SCPV) |l[ and is the subject of this article. 

The fact that spontaneous CP violation can occur in multi-Higgs extensions of the standard model is very well 
known 0. Typically, it is highly dependent on the structure of the Higgs potential. For instance in the Standard 
Model (SM) which has only one Higgs doublet, there is no physical phase in its vacuum expectation value (VEV) and 
hence no spontaneous CP violation. With two Higgs doublets [l[ , the ground state can break CP only if there are 
both quartic {(j)\(j)2) 2 as well as quadratic terms (p[<j>2 present in the potential. When cither of them is absent, the 
C ■ ground state corresponds to real VEVs of the Higgs doublets and there is no spontaneous CP breaking. In three Higgs 
extensions of theSM, even if the quadratic terms are absent, there can be spontaneous CP breaking0,|3| — manifested 
through both the exchange of charged [|[ or neutral [f| Higgs bosons. Since there is a widespread belief that the new 
scale physics includes supersymmetry, it is of interest to study CP breaking in MSSM and its multi-Higgs extensions. 

It is well known that in supersymmetric theories, holomorphy of the superpotential, which is a requirement for the 
theory to be supersymmetric, considerably restricts the terms in the Higgs potential. It is therefore important to 
study whether in such theories, spontaneous CP violation can occur. Such studies have indeed been carried out in 
simple extensions of MSSM and it has been shown that if we do not allow fine-tuning relations between the parameters 
of the superpotential, it is impossible to have spontaneous CP violation with only one additional singlet added to 
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MSSM_6] or with four and six Higgs doublet superfields (two or three pairs of H u ,Hd) extensions of MSSM 7]. The 
four-Higgs case remains incapable of spontaneously breaking CP even if one singlet is added provided that there are 
no dimensionful parameters in the superpotential involving the singlet Q; otherwise SCPV is in fact possible Q. 

The proof in the four- and six-Higgs case given in 0] uses geometric constructions. In this paper, we provide an 
alternative algebraic proof of the same result for the four Higgs doublet extension of MSSM and further show that 
our technique has the advantage that it can be extended to the N-Higgs extension for arbitrary N, with the result 
that there is no spontaneous CP breaking even in this case. We then show that this result remains true even in the 
presence of spontaneous R-parity breaking, where the sneutrino fields acquire nonzero VEVs[9(. 
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This paper is organized as follows: in Sec. [H] we discuss the general strategy for treating this question. In Sec. lIIIl 
we apply this to the general N-HMSSM, which has N/2 pairs of (H u , H d ) fields, and to the general N-HMSSM with 
spontaneous R-parity breaking; in particular, the 4-HMSSM is treated in more detail in Sec. lIII Cl as an illustration 
of the method. We present our conclusions in Sec. lIVI 



II. GENERAL MULTI-HIGGS MODELS AND CONDITIONS FOR SPONTANEOUS CP VIOLATION 

We consider extensions of standard model with N Higgs doublets <j> a with a = 1, . . . , N, denoted here as N-Higgs- 
doublet- models (NHDMs). In general we can assume that only the neutral members of the Higgs doublets acquire 
VEVs, which we parameterize as: 

(0 a ) = (O,^) T . (1) 

We can fix one of the phases 8i = from hypercharge invariance. Then, there will be SCPV (in the real basis) if 
5 a 7^ 0, 7r for some a > 2 and nonzero v a . Such violation breaks the canonical CP transformation 

a (i,x) C 4 P C(*,-x). (2) 

For some cases, there might be other inequivalent CP transformations that could remain as a symmetry (multiple CP 
symmetries). We do not treat this case here. 

We are interested to study the minima of the potential involving the N Higgs doublets and their CP properties. To 
build a SCPV model, it is necessary to begin with a CP invariant potential before SSB. There is a basis where CP 
invariance implies that all parameters of the potential are real when written in terms of (fy a 4>b- This is the real basis. 
Let us suppose that a CP invariant NHDM potential in this basis has the form 

V = V 1 (\^ a \ 2 ) + |0 a | 2 = ^ a . (3) 

Hence, we are avoiding terms such as |<?!>i| 2 0{</>2. We also assume we are seeking neutral minima and terms such as 
|(/>{</>2| 2 are equivalent to |</>i| 2 |</>2| 2 , hence, we include them in V\. 

This form of the potential covers a large class of theories such as the Weinberg model j^L H| for SCPV and supersym- 
metric models where the terms (fJ a 4>b, a ^ b, are only contained in the quadratic parts [7|. This can be also adapted 
for potentials, without trilinear terms, made up of the neutral components of different multiplets. 

Let us derive the extremum equations for Q, 

dV , BVi , dW n 

k^T^+> ,4>bk^T7- =0. (4) 



d<j>: k ^ aK dK aa ^j 0K dK ba 



We have used the shorthand 



If we contract ((U) with <jf ak , we obtain 



tk 



k=l,2 TaK b^a 

The imaginary part of ([6]) yields 



K ba = ^ h . (5) 



dV dV x „ dW 



dW 



= 0. (7) 



Spontaneous CP violation (SCPV) will be possible only if a solution for (JJJ is nontrivial, i.e., 

lm(K a ii) 7^ for some a ^ b, (8) 

without the need for any fine-tuning of parameters. 

To comply with a successful EWSB, it is necessary that 



K 246Ge V 



^K aa n — , (9) 
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which requires at least one of K aa non-null. Let us choose, without loss of generality, K\\ ^ 0. 

In the following, we will restrict the part of the potential that depends on the relative phases of the VEVs, i.e., W, 
to be either quadratic or quartic in the doublet fields. The quadratic case includes all general terms but the quartic 
case excludes further terms such as <j)\ 020103- To specify further the class of NHDM potentials we are considering 
here, we use the following parametrization for the potential: 



V = J2^a\ 2 + \ 5> Qb |0 a | 2 |0 b | 2 + W . 

a a, b 

The parameters (matrix) c ab are real and symmetric. 

A. Quadratic W 



(10) 



Let us choose 



W = ^ A Qb0l0b + h.C. = ^ ^abK ba 



(11) 



a<b 



a , b 



The parameters X a b (matrix) are real and symmetric, X ba = X a b, but we adopt X aa — 0. The terms in pT|) arc 
generic quadratic terms containing <jy a 4>b, a ^ b. We will refer to this choice as quadratic W . It is the case of some 
supersymmetric models 7]. For special potentials with 4, 5 and 6 Higgs doublets, it was shown that SCPV is not 
possible 0- We are interested in generalizing this result to an arbitrary number of doublets. 
First, we compute the derivatives for (fTTj) : 



dW 
dK, 



Xab 



bo 



Then, Eq. © yields 



Its imaginary part gives 



K„ 



dVi 
dKTa 



^ A afc -ftf ba = , a = 1, 



,N. 



^X a b^(K ba ) = 0, = 1, 



Instead of considering the real and imaginary parts independently, we can rewrite (|13[) as 

2 dVx 

b 

where we have used 



J2 X ^e Wba v a Vb = 0, a=l,...,N, 



K a b = v a v h e l6ah , v a = y/K aa , 9 ab 
After some straightforward manipulation we can rewrite (I15[) as 

dVi 

,e ~"\v a v b 

b 



J2e-^(v a v b -^S ab + X ab v a v b y Wb =0, a = l,...,N. 



(12) 



(13) 



(14) 



(15) 



(16) 



(17) 



In matricial notation the equation is equivalent to 



( 1 ^ 




(°\ 








W 



(18) 
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where the N x N real symmetric matrix B is defined by 



)ab = V a V b i^f 1 S ab + XabVaVb ■ 
OK na 



(19) 



Explicitly, 



■ ■ ^l,NVlV N \ 



V 



9Vi 



(20) 



NdK NN 



Equation (|18p implies that B as a function of {v a } should be a singular matrix. 

A noteworthy consequence of (I18p should be pointed out: the N equations define TV polygons in the complex plane. 
Each polygon, corresponding to the a-th row of (TT51) . is made up of N sides |(B) a i|, |(B) a 2|, • ■ • , |(B) a jv|, with (external) 
angles 6i,i-i — On — between the i-th side |(B) t| and the extension of the |(B) a ^_x| side; if some of (B) a j is 

not positive, we need to consider the angle On + tt instead of 9n. Thus we have a nontrivial solution (SCPV) if we 
can find a set of N polygons with the same angles Notice that among of the ^N(N — 1) sides (B) a j, only N 

of them can vary independently due to the independent {v a }- 



B. Quartic W 



Let us now choose 



(21) 



a.b 



a , b 



where Ab a = \ a b, & 7^ b, is real and X aa = 0. We will refer to this choice as quartic W. We assume terms such as 
|(^(/>k| 2 = |i^ a b| 2 are equivalent to K aa K bb for neutral VEVs and we include them into V\. Notice we are excluding 
terms such as <fi\4>2<j>\(f>3 or <j)\(f>2<p\4>i- This can be achieved by imposing a (Z-i) N symmetry on the potential. In this 
case 

dW 



Equation ^ is now 



dK, 



XnhKl 



bJ^ba ■ 



ba 



Kn 



dVi_ 



Ys^abKl^O, a = l,...,N. 



(22) 



(23) 



If we use (Unj) , we obtain 



a \dK aa 



0. 



(24) 



After some manipulations we can read 



(25) 



We can rewrite the last equation as 



( 1 \ 



1 » 



w 



(26) 
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where 



BV 

(MU = v 2 a —^5 ab + \ ab vy b . (27) 

In complete analogy with Eq. (fT8)) . Eq. ((26|) defines JV polygons with the only difference that now the angles are 
20 a i and the off-diagonal terms M ab depend quadratically on the VEVs. The JV sides of the a-th polygon are 
|(B) a i|, |(B) a|, ■ • ■ , |(B) a jv|, with angles 2# ij ,_i between the i-th side |(B) ai | and the extension of the |(B) 0) i_i| side. 
We may have a nontrivial solution (SCPV) if we can find a set of N polygons with the same angles 28 a \ and sides 
defined by (B) Qi . 

C. Constraints from nontrivial angles 



Let us analyze the polygon equations (|18j) . and, equivalently, f|26|) . The equations themselves imply B must be 
singular: 

rankB < JV - 1 . (28) 

However, if at least one polygon angle is nontrivial (#21 7^ 0, w) the singularity is such that 

rankB < JV — 2 . (29) 

In fact, one nontrivial angle implies at least another nontrivial angle. The reason for (|29[) is that we can write the 
first column of B as a linear combination of the N — 1 other columns as consequence of the real part of ([18]) . For the 
same reason, the imaginary part of (|18p — which does not depend on the first column of B — allows us to write the 
second column of B in terms of the N — 2 other columns to the right of the second column. 

The equality in Eq. (|28[) can be recovered in the absence of SCPV where all e l9il = ±1, i.e., a degenerate polygon 
confined to the real line. In this case, there is no relation among the rows of B, except from (|28p . Therefore, the 
JV equations corresponding to the JV rows of B are independent and can be used to find the JV quantities v a . This 
fact summarizes the relation between nontrivial angles (CP phases) and constraints on the sides (VEVs). These 
relations constrain the moduli v a and ultimately determines them. If there is not enough structure in the potential to 
allow nontrivial polygons, then there would be no CP breaking extremum. In particular, for N = 2 with cither only 
quadratic or quartic potential, a two-sided polygon always resides in the real line and then SCPV is not possible. 

For the simplest nontrivial polygon, i.e., a triangle, (|29[) is enough to find the nontrivial solution. To form a 
nontrivial triangle (finite area), we need rankB = 1. The geometrical counterpart of such a requirement is that 
triangles with the same external (or internal) angles must be similar. Conversely, similar triangles imply rankB = 1. 
This requirement can be translated to the fact that any 2x2 minor of the 3x3 matrix B is zero. If we take the 
relations coming from the vanishing of the three minors that only depend on one diagonal entry, then we get three 
equations that define V\ , V2, 1)3. Once we get a solution for the VEVs (sides of the triangles) we easily find the solutions 
for the CP phases (angles of the triangles). This solution is unique when it exists. 

Polygons with the same external angles are no longer guaranteed to be similar for quadrilaterals or polygons with 
more sides. The converse, however, remains true: if JV similar and noncontractible JV-gons are solutions for (|18[) . then 

rankB = 1 . (30) 

This case is only possible if all X ab have the same sign. The latter is just a necessary condition and, possibly for 
N > 3, a solution satisfying (1501) would require a special choice for the parameters of the potential; see Sec. lIIEl 

D. Algebraic consequences 

We have seen in (f2"9"]l that nontrivial CP phases require rankB < N — 2. Let us analyze its consequences. 
Let us assume rankB = N — 2. That means that N — 2 rows of B are linearly independent. Let them be the first 
JV — 2 rows. Since B is symmetric, we can take the following set as the (JV — 2)(JV + 3)/2 independent elements of B: 

M ai , a = 1, . . . , JV - 2, i = 1, . . . , JV (i > a) . (31) 

The elements B a b, with a, b > N — 2, can be written in terms of (|3Tj) as [10] 

M ba = X T b (M^-^y 1 Xa , (32) 



6 



where B^" 2 ) is the upper-left (TV - 2) x (TV - 2) submatrix of B and 

X a = (B la ,B 2a ,--- ,Bw-2,a) T , a,b=N-l,N. (33) 

Equations (|32jl and (|33[) are the key equations in our discussion of whether a multi-Higgs supcrsymmetric theory can 
support spontaneous CP violation or not. If rankB < TV — 2, more elements can be written in terms of a smaller set 
of independent elements. Equations (|3~Tj) and (|3"2"|) change accordingly. The extreme case of rankB = 1 is trivially 
satisfied by (1521) . 



E. Avoiding fine-tuning 



Let us show here some situations that could lead to a fine-tuning of the parameters of the potential unless they are 
imposed as a consequence of symmetries. For that, we should notice that the off-diagonal elements of B depend on 
the VEVs either as \ a bV a Vb (quadratic W) or X a bVa v b (quartic W). The diagonal elements of B depends on the VEVs 
in more complex ways. Let us assume quadratic W to be specific; we only need to replace v a — > v\ in the terms that 
contain X a b- 

We know that for triangles (TV = 3) Eq. (f2"9")l implies there is only one possible solution for Eq. (|32l) . Let us think 
in the case of the next simplest case: TV = 4. To have polygons with nontrivial angles, we need rankB to be 1 or 2. 
If rankB = 1 — the case of similar quadrilaterals — all rows of B are proportional to the first row. Then any minor 
of order two (2x2 submatrix) is null. In particular, we can find a 2 x 2 submatrix above the diagonal that does not 
contain elements from the diagonal of B itself. Such a minor is 

det (££:)=*«(SA;;)"™"*= - < m » 

Therefore, if we require that all the VEVs v a are nonzero, we are forced to have 

det(^ 13 * 23 )=0. (35) 

\/\l4 A 2 4 J 

Such a relation can be satisfied only by fine-tuning the parameters unless one finds some symmetry which guarantee 
this. Barring such possibilities, one can conclude that the only possibility for SCPV to occur for the case of four Higgs 
doublets is to have rankB = 2; it is the maximal value where SCPV is possible. We still have the possibility to set 
some VEVs to zero. But that means that the problem reduces effectively to considering an 3 x 3 submatrix of B, i.e., 
triangles; see Sec. Ill Fl 

We can easily extend (f3"4")) to an arbitrary number TV of doublets. Many submatrices of the type of (|3"4")) can be 
constructed by drawing a square inside the upper-right triangle of off-diagonal terms. We can see that the largest of 
such submatrices has size n x n for TV = 2n or 2n + 1. Therefore, 

no fine-tuning and nonzero VEVs =>■ rank B > n , for N — 2n or N — 2n + 1 . (36) 

This condition, together with (|2T)|) . reduces the number of cases we have to treat. In particular, for TV > 3, similar 
polygons can not be solutions to the extremization problem (TT5)) [and (f2"6"|) ] if we avoid fine-tuning or null VEVs. 



F. Null VEVs 



Let us treat the case of null VEVs. Suppose vn — 0. Then the matrix B reduces from a TV x TV matrix to an 
effective (TV — 1) x (TV — 1) matrix B' because the TV-th row and the TV-th column of B is now entirely null. If the rest 
of the VEVs v a are nonzero, we can use the argument used to obtain the constraint ([2"§f and apply it to the effective 
non-null submatrix B'. We conclude that 

rankB' < (TV- 1) -2; (37) 

the rank of B' is reduced by at least two units and the same conditions of Sees. Ill Dl and UlEl apply as if the relevant 
quantity were B'. If a number p of VEVs are null, we can extract the effective submatrix B' of size m x m, m — TV — p, 
by eliminating the p rows and columns made entirely of zeros. If we require at least one nontrivial relative angle 
(corresponding to two nonzero VEVs), we still need 

rankB = rankB' < m - 2 . (38) 

The condition (|36|) also changes to 

no fine-tuning ==>■ rankB > n , for m = 2n or m = 2n + 1 . (39) 
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III. APPLICATIONS TO NHMSSM 



Let us show here that there can be no SCPV at tree level in N-Higgs-doublets supersymmetric extensions of the 
SM if (i) we avoid fine-tuning of the parameters of the potential and (ii) the quartic terms come solely from D-terms. 
The proof is essentially an application of the constraints explained in Sec. Ill Dl 

We consider two cases: (A) an arbitrary number of N = 2n Higgs doublets and (B) rif = 3 effective Higgs doublets 
coming from the sleptons Li, i = 1, 2, 3, in addition to an arbitrary number of 2n Higgs doublets; the effective number 
of doublets is N — 2n + rif. The latter case is usually referred to as R-parity violating models [7]. 

In both types of models, we label the n + rif (n/ = or n/ = 3) doublets with hypercharge Y = — 1 with indices 
a = 1,.. .,7i + 7j/, and the n doublets with hypercharge Y = 1 with the indices a = n + rif + 1, . . . , N; N — 2ri + rif. 
Moreover, for Y = 1 doublets we just denote H a — ¥ (f> a while for Y = — 1 doublets we associate H a — > (f) a — i<Ji4>* a - All 
the effective doublets now have hypercharge Y = 1 and they are all denoted by {4> a }, a = l,...,N;N = 2n + rif. 



A. N = 2n Higgs doublets 

Let us consider N — 2n Higgs doublets. We are assuming that the quartic part of the potential comes only from 
the D-term 

V D = \G[\^\ 2 + |<fe| 2 + ■ ■ ■ + |0„| 2 - |0„ +1 | 2 |0 2 „| 2 ] 2 , (40) 

where G = {g 2 + g' 2 )/^. We also assume the VEVs of the doublets do not break electric charge. The potential (|4H|) 
is translated, in the notation of (fTO)) . into 

(Cab) = G ~tf) ' ( 41 ) 

where A n = A nxn is an n x n matrix with (A n ) a i, = 1 for all entries. The general potential will have the form (|10[) 
with quadratic W as (flTI) . Except for (c a t), which is defined in (14"T1) . all parameters of the potential are generic as 
long as they can provide a neutral vacuum and a potential bounded from below. 

One peculiarity of real potentials with quadratic W is that we can rewrite B in Eq. (|T9]) as 



l = diag(«i,i;2 ! ---,UJv)B diag(«i, v 2 , . . . , v N ) , (42) 

where 

B g fr = Sgbifig + ^ C ad V%) + Aqfe . (43) 

d 

The matrix B depends on the VEVs v a only in the diagonal terms and such a dependence occurs only through the 
terms ^^CadV 2 ,. In our particular case, the matrix (c a b) in Eq. (|4ip is so restrictive that it allows us to write 

Y,ca d v 2 d = e a f(v 2 ), (44) 

d 

where e a = 1 for 1 < a < n and e a = — 1 for n < a < N . The quantity f(v 2 ) is the combination 

n 

f{v 2 )^G^l~vl +a )- (45) 

o=l 

We will see that the dependence of B only on the combination f(v 2 ) is too restrictive to allow SCPV, unless a 
fine-tuning of the parameters is allowed. 

We first allow the possibility that some number p > of VEVs are null, and no more. Let us eliminate the null 
rows and columns and extract the effective submatrix B' of size m x m, where m = N — p. The requirement of at 
least one nontrivial CP phase imposes (|3"5)) . i.e., 

r = rank B = rank B' <m-2. (46) 

The possibility that B' has more null rows (columns) is excluded because that would require either more null VEVs 
or entire rows of (X a b) to be null. 
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Let us rearrange the rows (columns) of B in such a way that the first m rows (columns) correspond to the rows 
(columns) of B'. We also rearrange the rows (columns) of B' in such a way that the first m — 2 rows (columns) are 
linearly independent. We relabel (i a , c a b, X a b accordingly. All this rearrangement maintains the diagonal entries of B 
in the diagonal and we still have 



%aa = vl[^a + £af(v 2 )] , a = 1, . . . , 171 , 



(47) 



where f(v 2 ) is given in (|4"5")) and e a = ±1. 

Let us suppose now that m > 3. Nontrivial CP phases are not possible for m = 1 or m = 2, even for generic 
potentials with quadratic or quartic W. Take B' r ' as the r x r upper-left submatrix of B (or B'); r < m — 2 by (14^1) . 
The matrix B^ r ^ is nonsingular by construction. Then a generalization of Eq. (|32[) applies at least to B mm , B m _i. m _i 
and B TO _x m , i.e., 



where 



We can also define \ a as 



so that 



$m-l,m-l — xI-l( B ) 1 Xm-l 

® m>m =x r m (® {r) r l xm, 

B m -l,m = Xm-l( B(r) )" 1 Xm, 



Xa = (Bi a ,B 2a , ■ ■ • ,B r!a ) T , a — in 1 or vn . 

Xa = (Ala, A2a, ■ • ■ j ^r,a) T j 

Xa = w a diag(ui,« 2 , ■ ■ .,«r)Xo , a = m-l,m. 



Then Eq. pal) can be rewritten as 



t4_lB m _i jm _i - Xm-l( 



} ) Xm-1 : 



^X™(B (r) ) _1 X™, 



(48) 

(49) 
(50) 
(51) 

(52) 



where 



4W = 



A21 ^i+e 2 f{v 2 ) 



Air 
A2r 



(53) 



\ A n i A r , 2 ••• Hr + e r f(v 2 )J 

Now we use the fact that u TO and v m -i are non-null by hypothesis. We then obtain 



Mm-l + £m-l/(w 2 ) = Xm-l(® (r) ) Xm-1 , 

r 

xl.^y'xr 



(54a) 
(54b) 
(54c) 



Notice that these three equations only depend on the VEVs through the combination in f(v 2 ). Therefore if we solve 
Ea. (|54ap for f(v 2 ) in terms of the parameters of the potential, Eqs. (|54b[) and (|54cj) will give us two relations between 
the parameters of the potential. Such relations will be only satisfied by fine-tuning the parameters and is not stable 
under radiative conditions. Therefore, we conclude that CP phases must vanish and hence there is no spontaneous 
CP violation. 

If r < rn — 2 strictly, we can add to the list of equations in (|5"4"]) m — 2 — r more relations coming from Eq. (|3"2")l 
applied to the diagonal entries B aa , r < a < m — 2. We can also add more relations associated to the off-diagonal 
entries B a t,, r < a < m — 2, r < b < m— 2, a < b. Then more relations have to be satisfied simultaneously by the only 
combination in f(y 2 ) and the fine-tuning is worsened. 

Therefore, we conclude that the quartic terms — coming from D-tcrms — in supersymmetric multi-Higgs-doublet 
extensions of the SM do not have enough structure to break CP spontaneously, independently of the number of 
doublets N = 2n. 
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B. R-parity-violating models 

Let us consider N = 2n + rif doublets, rif — 3 coming from the sleptons Li, one for each family i; see Ref. 
We can still assume that the quartic terms, for neutral VEVs, comes solely from the D-term 



v d = \g[ E i^i 2 - E i^i 2 ] ■ ( 55 ) 

a—1 a—n-\-n f 

The matrix (c a b) is now 

(c ah ) =C(_ - A (n+ n/ )xn\ (56) 

As in (|41l) . the matrices A nxm are n x m matrices where all entries are unity. 

Now, the diagonal entries of B in (l4"3"f still depend only on one combination of the VEVs, 

^2c ad vj = e a f(v 2 ) , (57) 

d 

where now 

n+rif N 

/(, 2 )^g[e^- e A- (58) 

a— 1 a— ra+n/ + 1 

The factors c a are 6 a — -1-1 for 1 *C o, *C 71 -\- Ttf and 6 a — 1 for tl -\- Ttf <C cl A. Therefore the arguments used 
in Eqs. (f4"6")) - (|54cj) still apply and it is not possible to find solutions to the extremum equations with nontrivial CP 
phases. Thus it is not possible to have SCPV in supersymmetric multi-Higgs extensions of the SM with spontaneous 
R-parity violation. 



C. 4-HMSSM 

In this subsection, we illustrate the application of our technique to the case of 4-HMSSM (i.e., 4-Higgs-doublet 
supersymmetric extension of the SM) [7j . In this case, the quartic part of the potential (|40j) is simply 

Vb = iG[|^| 2 + |0 2 | 2 -|0 3 | 2 -|^| 2 ] 2 . (59) 
This is equivalent to considering the matrix (c a b) as 

/ 1 1-1 -1\ 

(60) 

-11 1 / 

where G = (gf + gf)/4. The function / (|4"5"]l in this case is 

f{ v *)^G[v{+vl-vi-vi]. (61) 

Allowing for null VEVs, we need at least three non-null VEVs to have at least two nontrivial CP phases. Therefore, 
the effective non-null matrix B' would have size m = 3 or m — 4. For both cases, nontrivial CP phases require 
r = rankB < m — 2. That in turn, implies the relations (|5l)) . 

For m — 4 (no null VEVs) and rankB = 2, considering the first two rows of B to be linearly independent, we have 
explicitly 




ft 2\ (\ \ \ ft Jt l + f( v2 ) ^12 \ 1 

M4 -/(«) = (A M A 24 ) ^ Xn M2 + /(w2) J 



- 1 / A 



13 
A23 

A14 
A24 



(62) 



A 34 -(A 13 A 23 )^ ^f^) [x 24 
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For Tji — 4 and rankB = 1, (|rJ2"j) is replaced by 



» + ^ = lZ+fWV (63a) 

- f(v 2 ) = — ^frrr • ( 63b ) 



>i + /(« 2 ) 

~Ml + /(« 2 ) 



M4 - = ^7 . ( 63c ) 



and 



\ A12A13 

A23 = , tl 2 , , 64a) 

Mi + /(« 2 ) 

A24 = i 77 2\ ' 64b ) 

Ml + /(^) 

A13A14 

^34 = — - 77 2 ^ ■ (64c) 
Mi + /(« ) 



For m = 3 and rankB = 1, assuming = 0, the relations for B22 and B33 are the same as Eqs. (I63a[) and (163b|) ; 
the relation for B23 is the same as (|64a[) . If other VEVs are null instead of V4, we obtain similar relations. 
In all cases, fine-tuning is necessary to satisfy the requirements of nontrivial CP phases. 



IV. CONCLUSIONS 



To summarize, we have seen that the possibility of SCPV in NHDM potentials where complex combinations <j>\<f>b 
appear only on either quadratic (<j>' a (f>b) or quartic (4>t4'b) 2 parts can be related to the existence of N nontrivial polygons 
defined by a N x N real matrix B which is a function of the moduli \<fi a \ — v a as in (flT)|) or (f2~T|). These polygons 
generalize the triangle described in [3| for the Weinberg 3-Higgs-doublet model. For models with a moderate number 
of doublets, the characterization of the possible extrema through polygons should be still useful to guide the search 
for solutions. 

As a direct application of these ideas, we have shown that in arbitrary Higgs-doublets extension of MSSM, it is 
impossible to break CP spontaneously at tree level unless we allow fine-tuned relations between tree level parameters 
of the theory. This result remains true even if the sneutrinos acquire VEVs so that R-parity is spontaneously broken. 
This generalizes the earlier result in [7[ for two or three pairs of Higgs doublets to the case with an arbitrary number of 
them. A simple corollary of this result is that in supersymmetric left-right models with arbitrary number of bidoublets 
where B — L is broken by B — L = 2 triplets, there is no spontaneous CP violation in the bidoublet sector. This result 
has application to a solution of the strong CP problem without the need for an axion (see, for instance, [ill]). 
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